arXiv: 1509.00213vl [hep-lat] 1 Sep 2015 


Renormalization functions for Nf=2 and 


Nf=4: Twisted Mass fermions 





C. Alexandrou“’^, M. Constantinou“’*', 


H. Panagopoulos“El 


a Department of Physics, University of Cyprus, 
POB 20537, 1678 Nicosia, Cyprus 


b Computation based Science and Technology Research Center, 
The Cyprus Institute, 

15 Kypranoros Str., 1645 Nicosia, Cyprus 


We present results on the renormalization functions of the quark field and fermion bilinears with 
up to one covariant derivative. For the fermion part of the action we employ the twisted mass 
formulation with Nf=2 and A^/=4 degenerate dynamical quarks, while in the gluon sector we use 
the Iwasaki improved action. The simulations for A^/=4 have been performed for pion masses in 
the range of 390 MeV - 760 MeV and at three values of the lattice spacing, a, corresponding to 
/3=1.90, 1.95, 2.10. The Nf=2 action includes a clover term with Caw=l.57551 at /3=2.10, and three 
ensembles at different values of m-,,. 

The evaluation of the renormalization functions is carried out in the RF scheme using a momentum 
source. The non-perturbartive evaluation is complemented with a perturbative computation, which 
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I. INTRODUCTION 

Over the last years, simulations of Quantum Chromodynamics (QCD) have advanced remarkably and are, 
nowadays, being carried out at close-to-physical values for the parameters of the theory. Therefore, ab initio 
calculations of hadron structure within lattice QCD yield results that can be connected to experiment more 
reliably than ever before. A number of lattice groups are producing results on nucleon form factors and first 
moments of structure functions at or close to the physical regime both in terms of pion mass as well as in terms 
of the continuum limit (see Ref. [l| and references therein). At the same time properties of other hadrons that 
are difficult to study experimentally are being pursued within lattice QCD. These include the axial charges 
of resonances such as the A Q or other nucleon excited states Q, hyperons [1-Q and charm baryons [y]. 
For all these quantities, one needs the renormalization functions in order to obtain the continuum predictions. 
Moments of generalized parton distributions (GPDs) are connected to generalized form factors and provide 
detailed information on the internal structure of hadrons in terms of both the longitudinal momentum fraction 
and the total momentum transfer squared. Beyond the information that the form factors yield, such as size, 
magnetization and shape, GPDs encode additional information, relevant for experimental investigations, such as 
the decomposition of the total hadron spin into angular momentum and spin carried by quarks and gluons. In 
lattice QGD one calculates matrix elements of fermion operators between the relevant hadron states and unless 
these operators correspond to a conserved current they must be renormalized in order to extract the physical 
information one is after. In many cases, calculation of renormalization functions (RFs) can be carried out using 
lattice perturbation theory, which proves to be extremely helpful in cases where there is a mixing with operators 
of equal or lower dimension, such as the chromomagnetic operator [3, Q and the operator measuring the glue of 
the nucleon ilS- However, perturbation theory is reliable for a limited range of values of the coupling constant, 
g, and of the renormalization scale, g. For this reason, a non-perturbative computation on the RFs is preferable. 

In this work we will combine both perturbative and non-perturbative computations in order to obtain an 
improved evaluation for the RFs of the quark field, ultra-local and one-derivative fermion operators within the 
twisted mass formulation of Wilson lattice QGD [^. In particular, we compute lattice artifacts to all orders in 
the lattice spacing, a, using one-loop perturbation theory and we subtract them from the non-perturbative results 
for the RFs. This subtraction suppresses lattice artifacts considerably depending on the operator under study 
and leads to a more accurate determination of the renormalization functions. We show that lattice artifacts are 
non-negligible in most cases, and are significantly larger than statistical errors. 

We use the Rome-Southampton method (RF scheme) [l^ to compute the renormalization coefficients of arbi¬ 
trary quark-antiquark operators non-perturbatively. In this approach the renormalization conditions are defined 
similarly in perturbative and non-perturbative calculations. The RFs are obtained for different values of the 
renormalization scale, and on several ensembles corresponding to different pion masses, so that the chiral limit 
can be safely taken. Since the goal is to make contact with phenomenological and experimental studies, which 
almost exclusively refer to operators renormalized in the MS scheme, one needs the renormalization functions 
leading from the bare operators on the lattice to the MS operators in the continuum. The conversion to the MS 
and the evolution to a reference scale of 2 GeV is performed using three-loop perturbation theory. 

The paper is organized as follows: Section [IT] presents the lattice formulation and gives details on the gauge 
configurations and the parameters of each ensemble. Section |IT] includes the definition of the operators under 
study, as well as the renormalization conditions for the RT scheme. The methodolody of the non-perturbative 
computation is described in Section IIIII Section focuses on the perturbative procedure for the evaluation of 
the one-loop lattice artifacts to all orders in the lattice spacing denoted by 0{g‘^ a°°). This is a crucial component 
of this work, since the subtraction of the 0{g'^ a°°) contributions from the non-perturbative estimates of the RFs 
leads to removal of the bulk of lattice artifacts. The main part of the paper is Section El which presents the 
results of this work, including their chiral extrapolation, the conversion to the MS scheme and the evolution 
to a reference scale of 2 GeV, via the intermediate Renormalization Group Invariant scheme. Particular focus 
is given to the a°“)-corrected data for the RFs, and we show, for selected cases, a comparison with the 
0(5^ a^)-corrected expressions. The final values of the chirally extrapolated RFs at the limit (ap)^ —>■ 0 are 
presented in Table Hill In Section IVIII we give our conclusions. For completeness we provide in Appendix lAl all 
necessary formulae for the conversion to the MS scheme, and in Appendix IbI we present the 0{g‘^ a°°)-corrected 
RFs for all the previous twisted mass fermions ensembles fl^. Il3|. recomputed in the framework of this work. 


II. FORMULATION 
A. Simulation details 

The gauge field configurations were generated by the European Twisted Mass Collaboration (ETMC) employing 
the twisted mass fermion action. We use ensembles generated with V/=4 light degenerate quarks [l^ with the 
twisted mass action, as well as Nf=2 degenerate quarks 0 , in which a clover term is also included in the fermion 
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action with Csw=l-57551. We note that the renormalization functions computed using the Nf=A ensembles will 
be applied to renormalize matrix elements computed using Nf=2+1+1 gauge field configurations [I3- We adopt 
the Rh renormalization scheme, which is mass independent, and consequently the RFs are defined at zero quark 
mass. For this reason, the A^/=2+l+l ensembles cannot be used to compute the RFs due to the fact that the 
mass of the strange and charm quarks are fixed to their physical values, and extrapolation to the chiral limit is 
not possible. Therefore, in order to compute the renormalization functions needed to obtain physical observables, 
ETMC has generated iV/=4 ensembles at the same values of /3 so that the chiral limit can be taken. Details on 
the simulations can be found in Refs, [l^ [l8j|. 

Automatic 0{a) improvement for twisted mass fermions may be achieved with maximal twist, by tuning the 
mpcAC quark mass to zero. For the case of A^/=4 configurations, achieving maximal twist is difficult particularly 
if the lattice spacing is not very fine, which is associated with a change in the slope of mpcAC with respect to 
1/(2 k) [I^. In order to tackle this issue, Monte Carlo simulations are performed in pairs not exactly at maximal 
twist but with opposite values of mpcAC- As proposed in Ref. [2^, by averaging the RFs computed on ensembles 
with opposite values of mpcAC, 0(a) improvement is achieved (see also Ref. [2l| and references therein). 

In the gluon sector we use, for all ensembles, the Iwasaki improved gauge action [2^, which includes besides 
the plaquette term also rectangular (1 x 2) Wilson loops 

*^5 = f E E - ReTr(C/i,(^;j} + b^J2{l- ReTr(C/if J} ] (1) 

with j3=2 Nc/bi = —0.331 and the (proper) normalization condition b n = 1 — 86i = 3.648. 

The simulation details, the parameters and the values of the pion mass (23| of each ensemble used in this work 
are given in Tables HI-HIl for the A^/=4 and A'f = 2 ensembles, respectively. The values of the lattice spacing have 
been determined using the nucleon mass [ij, UM ■ 


afi 


c^Mpcac 

aKdps 

lattice size 

P = 1.90, a = 0.0934 fm 

0.0080 

0.162689 

-h0.0275(4) 

0.280(1) 

24® 

X 48 


0.163476 

-0.0273(2) 

0.227(1) 



0.0080 

0.162876 

-h0.0398(l) 

0.279(2) 

24® 

X 48 


0.163206 

-0.0390(1) 

0.241(1) 



/3 = 1.95, 0 = 0.082 fm 

0.0020 

0.160524 

-h0.0363(l) 


24® 

X 48 


0.161585 

-0.0363(1) 




0.0085 

0.160826 

-h0.0191(2) 

0.277(2) 

24® 

X 48 


0.161229 

-0.0209(2) 

0.259(1) 



0.0180 

0.160826 

-h0.0163(2) 

0.317(1) 

24® 

X 48 


0.161229 

-0.0160(2) 

0.292(1) 



j3 = 2.10, a = 0.064 fm 

0.0030 

0.156042 

4-0.0042(1) 

0.127(2) 

CO 

to 

CO 

X 64 


0.156157 

-0.0040(1) 

0.129(3) 



0.0046 

0.156017 

-40.0056(1) 

0.150(2) 

32® 

X 64 


0.156209 

-0.0059(1) 

0.160(4) 



0.0064 

0.155983 

-40.0069(1) 

0.171(1) 

CO 

to 

CO 

X 64 


0.156250 

-0.0068(1) 

0.180(4) 



0.0078 

0.155949 

-40.0082(1) 

0.188(1) 

32® 

X 64 


0.156291 

-0.0082(1) 

0.191(3) 




TABLE I: Simulation details for the A/=4 ensembles of twisted mass fermions. The lattice spacing is determined using 
the nucleon mass computed with A/=2+l+l twisted mass configurations at the same values of /3. 


The number of configurations in each ensemble varies between 10 to 50 separated by 20-100 trajectories, 
depending on the ensemble. The small size of these ensembles, is more than sufficient for use of the momentum 
source method, which offers high statistical accuracy, easily below 0.5% even for 10 configurations (see Section lllll) . 
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ap, K aMps 

lattice size 

P = 2.10, a = 0.093 fm, c, 

= 1.57551 

0.0009 0.13729 0.0621(2) 

48® X 96 

0.0030 0.1373 0.110(4) 

24® X 48 

0.0060 0.1373 0.160(4) 

24® X 48 


TABLE II: Simulation details for the Nf=2 twisted mass ensembles with a clover term. The lattice spacing is determined 
using the nuleon mass computed using the same Nf=2 ensembles. 


In our computation, we mostly use “democratic momenta” in the spatial direction, such as: 


(ap) = 2TT 



1 Tlx Tlx Tlx \ 

2Lt'lZ'lZ'Ts) 


where Lt (Lg) is the temporal (spatial) extent of the lattice and Tit and tIx take values within the range: 


( 2 ) 


nte[2,20], n,e[l,10], (3) 

depending on the lattice size of each ensemble, so that they correspond to momentum up to (ap)^~7. To 
fill in some gaps between the momentum ranges we also include a few non-democratic momenta of the form 
(ut, Uj;, u-x ±1), which show similar behaviour with neighbouring democratic momenta. 


B. Definition of operators and renormalization prescription 


In this work we consider ultra-local fermion operators: 


Os = XT“X 

Op = X75t“x 

O^ = 

Oa =X757A.r“x 

Op = XCTpi^T°'X 


I 

I — 


0 = 1,2 

0 = 3 

o = 1,2 
0 = 3 


a = 1 

-'0757m'^V 0 = 2 
0 = 3 


0 = 1 

-^7/i'rV 0 = 2 

V'757m'^^V' a = 3 




o = 1,2 
0 = 3 


Opp = X750-^i-'r“x = 


and the following one-derivative fermion operators: 


= X7{/.^i.}T“X 

Oba^ =Xl5l{p'^i^}T°'X 
= X750-,x{i/^p}t“x 


V'75CTpi.T“V' 0=1, 
a = 3 

2 


D Ip 

0=1 

< -'075^p D 

0 = 2 

y’4’l{pD Ip 

0 = 3 


VX{pD^T^tp 

0=1 

< -'lpX{f,D ^yp^lp 
\ ^ ^ 

0 = 2 


0 = 3 

f— TT 

1 wiT>(yp{v i^}T°'tp 

o = 1,2 

\-iip(Tp{u D p^l-ip 

0 = 3 


(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 
(9) 


( 10 ) 


( 11 ) 


( 12 ) 
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all given in the 


where 

i 

and 

tlj{x)^ fj, = - 

a 

For completeness we include in the list even though its components are related to those of O^. We denote the 
corresponding RFs of the ultra-local fermion bilinears by Zg, Zp, Zy, Z'y, Z^, ^fp- ^ massless renormalization 
scheme, such as the RF, the RFs are defined in the chiral limit, where iso-spin symmetry is recovered. Hence, the 
renormalization functions become independent of the isospin index a = 1,2,3 and we drop the a index from here 
on. Still note that, for instance, the physical is renormalized with Zp^ while needs Zy, which 

differ from each other even in the chiral limit. 

The one-derivative operators are symmetrized over two Lorentz indices and are made traceless: 

Q{<yr} ^ (16) 

A 

which avoids mixing with lower dimension operators. The corresponding RFs of the one-derivative operators 
are denoted by Z^y, Z^p^, Z^j,. The one-derivative operators fall into different irreducible representations of 
the hypercubic group, depending on the choice of indices. Hence, we distinguish among them according to the 
following 


Odvi 

= Odv with p = V 

(17) 

C’dV2 

= Odv with p A ^ 

(18) 

C’dai 

= Oba with pL = V 

(19) 

C’dA2 

= Oda with p A ^ 

(20) 

OuTl 

= Obt with p ^ V = p 

(21) 


= Obt with p, ^ V ^ p ^ p . 

(22) 


Thus, .^DVi will be different from Zdv 2, but renormalized matrix elements of the two corresponding operators 
will be components of the same tensor in the continuum limit. 


sted and physical basis as shown above. The covariant derivative is defined as: 
- ^ 

2 I 2 2 


Uf_c{x)i’{x -I- a/t) - ipix) 


and '^*'ijj{x) = —- 


U'^Ax — ajl)'tp{x — ajj.) — 


(13) 

(14) 


i/jix + afi)Ul{x) -A{x) 


and i/;(x)V* = —- 


'ip{x — afl)UAx — ajl) — 


(15) 


The renormalization functions are computed in the RF scheme at different renormalization scales, The RFs 
are determined by imposing the following conditions: 


Zq = ^Tr[(5^(p))-i5«°“(p)] 


(23) 


^q-'^OY^Tr[r^(p)r 


Born-1 


(f)] 


= 1 . 


(24) 


where the momentum p is set to the renormalization scale p. The trace is taken over spin and color indices, 5'®°’’'’ 
is the tree-level result for the propagator, and F®”''*' is the tree-level expressions for the fermion operators, that 
is 


for the ultra-local bilinears, and 


O 


{p,v} 

DV 


1 

2 1 


F °™(p) = 1, yg, y^, ygy^, yg , 




(25) 

(26) 


O 


{nv} 

DA 


^'ysy^. D 


_ -O 


^ 757r 


D r'S 


c/dt 


^ ^ 75^^!. £> p 4' -k 4' ygCTpp D - i^^p ^ 4' 75 ^^^ D 


(27) 

(28) 
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for the one-derivative operators. The presence of 5'®°''" and ensure that Zq = 1, Zq = 1 when the 

gauge field is set to unity. The RT values for the RFs are translated to the MS scheme at ^ =2 GeV using an 
intermediate Renormalization Group Invariant scheme. 


III. NON-PERTURBATIVE CALCULATION 

For the non-perturbative evaluation we follow the same procedure as our previous work , and here we 

summarize the important steps of the calculation. We first write the operators in the form 

0{z) ='^u{z)J{z,z')d{z'), (29) 

z' 

where u and d denote quark fields in the physical basis and denotes the operator we are interested in. For 
example J'(z,z') = Sz,z'Jfj. corresponds to the local vector current. For each operator we define a bare vertex 
function given by 

= ^ X! e~"^^^~y'>{u(x)u{z)Jiz,z')diz')d{y)), (30) 

x,y,z,z' 

where p is a momentum allowed by the boundary conditions, V is the lattice volume, and the gauge average, 
denoted by the brackets, is performed over gauge-fixed configurations. The Dirac and color indices of G{p) are 
suppressed for simplicity. 

We employ the approach, introduced in Ref. [2^, which uses directly Eq. (l30|) without employing translation 
invariance and one must now use a source that is momentum dependent but can couple to any operator. For 
twisted mass fermions, we use the symmetry S'“(a:, y) = 'j 5 S‘^'^{y, between the u— and d—quark propagators, 
and therefore, with a single inversion one can extract the vertex function for a single momentum. The advantage of 
the momentum source approach is a high statistical accuracy and the evaluation of the vertex for any operator at 
no significant additional computational cost. The drawback is that we need a new inversion for each momentum. 
We fix to Landau gauge using a stochastic over-relaxation algorithm , converging to a gauge transformation 
which minimizes the functional 


F = ^ Re tr [U^[x) + U^x - p)] . (31) 

The propagator in momentum space, in the physical basis, is defined by 

S^p) = yJ2 , S^{p) = yJ2 {d{x)diy)) , (32) 

x,y x,y 

and an amputated vertex function is given by 

r{p) = {S-{p))-^G{p){S‘^{p))-^. (33) 

Finally, the corresponding renormalized quantities are assigned the values 

Sr{p) = Z^S{p), TR{p) = Z-^Zor{p) . (34) 

In the twisted basis at maximal twist, Eq. (l30l) takes the form 

^ iU ^ ((1 -b ij 5 )uix)u{z){l + ij5)Jiz, z'){t - yB)d{z') d{y){l - iqs)) , (35) 

x,y,z,z' 

which simplifies when using the relation between the u- and d-quark propagators, that is iS“(x, z) = '~f^S‘^\z, x)^^. 
After integration over the fermion fields it becomes 

G(,p) = X] “ ii5)S'^\z,p){t -yz)J{z,z'){t - i75)5'^(z',p)(l - fqs)^ , (36) 


^ In an alternative approach that relies on translation invariance, one may shift the coordinates of the correlators in Eq. (I30I I to 
position 2 = 0 fSHl . 
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where {■■■)^ denotes the integration over gluon fields, and S{z,p) = ^ye'^P^S{z,y) is the Fourier transformed 
propagator with respect to one of its arguments, on a particular gauge background. It can be obtained by 
inversion using the Fourier source 


bl{x) = e^^^5^p5ab, (37) 

for all Dirac a and color a indices. The propagators in the physical basis given in Eq. (|32ll can be obtained from 

Z 

(38) 

Z 

which only need 12 inversions (instead of 24) despite the occurrence of both u and d quarks in the original 
expression. We evaluate Eq. (1351) and Eq. (1381) for each momentum separately employing Fourier sources over a 
range of (ap)^ for which perturbative results can be trusted and finite a corrections are reasonably small. The 
amputated vertex functions of Eq. (13311 computed for each operator, as well as the inverse quark propagator, 
enter the renormalization prescription of Eqs. (USD - (I21- 


IV. ONE-LOOP CALCULATION OF ARTIFACTS TO ALL ORDERS IN THE LATTICE SPACING 

An improvement over previous work Bin, where we evaluated the 0{g'^ a^) perturbative artifacts, is the 
computation of the one-loop perturbative artifacts to all orders in the lattice spacing, 0{g^ a°°). These artifacts 
are unavoidably present in the non-perturbative vertex functions of the fermion propagator and fermion operators 
under study. In our previous work [mil, the 0{g'^ a^) perturbative artifacts were subtracted from the non- 
perturbative RFs, leading to improved estimates. However, for large values of the scale (ap)^, the 0{g'^ a?) terms 
tend to increase becoming, thus, unreliable. 

As will be demonstrated by our results, the lattice artifacts depend on the operator under study, as well as on 
several parameters such as the coupling constant, the fermion and gluon action parameters, the lattice size, the 
lattice spacing and the renormalization scale. Thus, a proper subtraction of the finite lattice size effects from 
the non-perturbative values requires a separate perturbative evaluation of the 0{g^ a°°) terms for each ensemble 
and each value of the four-momentum, in order to match our non-perturbative computation. In other words, the 
0{g^ a°°) contributions that are subtracted from each black circle point shown in the plots of Section IVl requires 
a separate perturbative computation. Unlike the case of the 0{g^ a^)-subtraction used in our previous work for 
the renormalization functions [Billlll, the 0{g‘^ a°°) contributions cannot be given in a closed form. 

There are six Feynman diagrams that enter the perturbative computation: Two for the fermion propagator and 
four for the fermion operators, as shown in Figs. [1]-[51 The operator insersion is represented by a cross. In this 
work we restrict ourselves to forward matrix elements (i.e. 2-point Green’s functions, zero momentum operator 
insertions). The Feynman diagrams are evaluated using our symbolic package in Mathematica, and details on 
the algebraic operations can be found in Ref. 
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FIG. 1: One-loop diagrams contributing to the fermion propagator. Wavy (solid) lines represent gluons (fermions). 


As a general strategy in our perturbative computations for the RFs we employ a variety of fermionic and gluonic 
actions, in order to obtain results that are applicable to simulations performed by various research groups. In 
this paper, however, we only present the results for the twisted mass action including a clover term. The latter 
is kept as a free parameter and can be sent to zero as required for the Nf=A ensembles presented in Table HI 
Although the clover parameter is treated as free throughout the perturbative computation, for our final estimates 
in the Nf=2 ensembles it is set to its tree-level value suggested by one-loop perturbation theory. Csw = 1. 

The computation of the 0{g'^ a°°) terms was first employed by the QCDSF Collaboration for Clover 

fermions and Wilson gluons, and was later generalized to include more complicated fermion and gluon actions 
[S^. The main difference between the computation of the 0{g^ a°°) and the 0{g‘^ a?') terms, is that the latter are 
extracted by performing a Taylor expansion with respect to a. The 0{g^ a°°) terms, however, cannot be given in 





8 


FIG. 2: One-loop diagrams contributing to the computation of the fermion operators. A wavy (solid) line represents 
gluons (fermions). A cross denotes an insertion of the operator under study. In the case of the ultra-local operators, only 
the upper right diagram contributes due to the absense of gluons in the vertices. 

a closed form in terms of a (since it is included in the propagators) and a separate calculation is performed for 
each value of the momentum, (ap)^. Of course, from the resulting expression one must omit the 0{g^ a^) terms 
since we are interested only in the lattice artifacts. Note that, in most cases, the latter contributions include 
logarithms log{a ))), which are also subtracted. In a nutshell, the lattice artifacts to all orders in the lattice 
spacing are computed in the procedure summarized in the following expressions, in which the 0{g^ a°) terms 
computed in Ref. [l^, are omitted: 


DZq{a,p) = {Vg{a,p) - Vq{0,p)) 


(39) 


where 


DZo{a,p) 


( Vg(a,p) 
\Vo{a,p) 


MhA] 

Vo{0,p)) 



Vg(«,p) = ^Tr [{S^{a,p)) 1 ^^““(p)] , 


Vo(a,p) = ^Tr[r^(a,p)rB°™-i(p)] 


(40) 


(41) 


and S^{a,p), T^{a,p) are the results up to one loop and to all orders in a. Finally, the perturbative 0{g^a°°)- 
terms are subtracted form the non-perturbative values 


Zf'--^p,a) = Zf{p,a) - DZgia.p) (42) 

Zf’‘^'^\p,a) = Z^^'{p,a) - DZo{a,p ). (43) 

In Figs. [3]-H] we plot DZq{a,p)/g‘^ and DZo{a,p)/g^ for the ultra-local bilinears corresponding to the Iwasaki 
improved action using a lattice size of 24^ x 48 and several values of (ap)^ within the range of 0-4. For comparison 
we also include the corresponding 0(p^a°)/p^ terms computed in Ref. [l^. Since a clover term is included in 
the Nf=2 ensembles, we consider both values: Csw = 0 (left figures) and Csw = 1 (right figures). An immediate 
observation is that momenta with the same (ap)^ lead to different lattice artifacts, which is expected, since 
beyond 0(a°) these terms depend not only on the length, but also on the direction of the four-vector p, due to 
the presence of Lorentz noninvariant structures, such as: 


2^4 _ EpaVp 
“ p2 - Y.pa^Pl ’ 


(44) 


appearing to 0{g^a?). It is also interesting to see that the lattice artifacts depend on the operator under study in a 
non-predictible way since in some cases {Zq, Zy, Zt) the inclusion of a clover term diminishes the artifacts, while 
in another (.^s) h enhances them. For the case of Zp^ and Zp the lattice artifacts for Csw = 0,1 are comparable. 
As expected, comparison between 0{g^ a^) and 0{g'^ a°°) for small values of the momenta ((ap)^ 1) reveals 

a very good agreement, since the 0{g‘^ a?') terms are the leading contributions of the lattice artifacts. For larger 
momenta the difference between 0{g‘^ a^) and 0{g'^ a°°) is more apparent, as will be discussed in Section IVl 
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FIG. 3: DZq{a,p) /as a function of (ap)'^ with the Iwasaki gluon action and a 24® x 48 lattice for Csw = 0 (left) and 
Csw = 1 (right). 



FIG. 4: DZs{a,p)/g^ as a function of {ap)'^. The notation is the same as for Fig. [S] 
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FIG. 5: DZp{a,p)/g^ as a function of (ap)^. The notation is the same as for Fig.O 
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FIG. 6: DZv{a,p)/g^ as a function of (ap)^. The notation is the same as for Fig. [3] 


0.05 

0.04 

0.03 

0.02 

0.01 

0.00 

- 0 . 01 , 


0.05 r 


0(g"a^) 




O(g^a') 

• 0(g'a‘"') 

. •. ■. “ 

0.04 


• 0(g^a‘"') 


! Mt* : i 


•.V 




«• I «,• • I ”• 

. *• .* I* t* ; * 


-.1- I : 


0 


(a p) 


0.03 

0.02 

0.01 

0.00 

- 0.01 


* ♦ • *« . • • ♦ •*?* 




'■I.-.. 





0 


(ap) 


FIG. 7: DZA{a,p)/g^ as a function of (ap)^. The notation is the same as for Fig. [S] 



FIG. 8: DZT{a,p)/g^ as a function of (ap)^. The notation is the same as for Fig. [S] 
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V. NON-PERTURBATIVE EVALUATION 
A. Chiral extrapolation 

In order to obtain the renormalization functions in the chiral limit we perform an extrapolation using a linear fit 
with respect to We find that the RFs obtained in this work have a very mild dependence on the pion mass for 
all ensembles. In fact, with the exception of a few small values of (ap)^, there is no visible pion mass dependence 
within the small statistical errors. Allowing a slope and performing a linear extrapolation with respect to 
the data yield a slope consistent with zero. Figs. 151 -1121 demonstrate the pion mass dependence of the RFs using 
the Nf=2 and A^/=4 ensembles at ,0=2.10. The statistical errors are too small to be visible. Figs. 1^ - [TT] provide 
a more general picture of the m^-dependence by displaying the RFs as a function of the renormalization scale 
(p2 _ while the two plots of Fig. IT^show the data at (ap)^ = 3 as a function of the twisted mass 
These plots show clearly that the slope of the fit is consistent with zero. 

In this discussion the renormalization function of the pseudoscalar density, Zp, has been excluded since there 
is pion pole contamination that needs to be taken into account. Thus, a polynomial fit with respect to the pion 
mass is not suitable. An approriate chiral extrapolation of Zp and the ratio Zp/Zg is discussed in the following 
Subsection. 



(a p)^ 

FIG. 9: Pion mass dependence of Zq and the RFs of the ultra-local bilinears for Nf=2 at 0=2.10 as a function of the 
renormalization scale. 



FIG. 10: Pion mass dependence of the RFs of the one-derivative vector and axial operators for Nf=2 at 0=2.10 as a 
function of the renormalization scale. 
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FIG. 11: Pion mass dependence of the RFs of the one-derivative tensor operator for Nf=2 at /3=2.10 as a function of the 
renormalization scale. 




FIG. 12: Pion mass dependence of Zq and the RFs of the ultra-local bilinears (left plot) and the one-derivative operators 
(right plot) for Nf—4 at /1=2.10 as a function of a at (ap)^ = 3. 


B. Pion-pole subtraction of Zp and Zp/Zs 


The correlation functions of the pseudoscalar operator have pion-pole contamination which needs to be treated 
carefully. In order to subtract the pole contribution we use a 2- and 3-paraineter Ansatz for the pseudoscalar 
amputated vertex function, Fp, of the form: 


= ap 


cp 

ml 


= 


1 2 Cp 

ap + bpm^-\ - 5 -. 

mi 


(45) 

(46) 


The fit parameters depend on both the momentum and the value of /3 (i.e. ap = ap{/3^p)), and thus we estimate 
them separately on each value of p and /3. We find that the coefficient bp is very small and competes with cp in 
the 3-parameter fit. In addition, they both carry large statistical errors, which result in a large error in the final 
determination of Zp once the term cp/m^ is subtracted from the pseudoscalar matrix elements: 


psub _ p Cp 

Ip —Ip- j. 

mi 


A way around this problem is to employ the 2-parameter fit of Eq. (1451) directly to the ratio: 

rp(p, m^) 


Vp(p, m^) = 


^q(p, mT,)C] 


RI', MS 


(p,2GeV) 


(47) 


(48) 


where Cp ’ {p, 2 GeV) is the conversion to the MS scheme and the evolution to a scale of 2 GeV. This fit allows 

us to obtain directly Z^® in the MS-scheme and at the chiral limit ^ from I/ap. In a similar manner, we obtain 
directly Zs/Zp from the coefficient ap computed from the 2-parameter fit of 


Vspip,mT,) = 


Tp{p, m^r) 


( 49 ) 


^ Alternative fit functions and their stability are discussed in Ref. ISOl . 
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As an example of the pion pole contamination and its subtraction we show, in Fig. 1131 Vp(p, m-jr) and Vsp{p, PiTr) 
using the Nf=2 ensembles at /?=2.10 for each value of the pion mass before and after the sutraction of the pion 

pole term, '^^2 ^ ■ Fig. [14] is similar to Fig. [T3| for the fV/=4 and /3=2.10 ensembles. The range of the y-axis 
is the same for the unsubtracted and subtracted cases in order to see clearly the effectiveness of the pion-pole 
subtraction. 
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FIG. 13: Upper panels: Eq. (I48II (left plot) and Eq. (I49II (right plot) for Nf=2 at /3=2.10 as a function of {ap)^ for each 
value of the pion mass before the pole subtraction. Lower panels: Eq. (I48II (left plot) and Eq. (I49II (right plot) for Nf=2 
at /3=2.10 as a function of (ap)^ for each value of the pion mass after the pole subtraction. 



FIG. 14: As Fig. 1131 but for the A^/=4 and /3=2.10 ensembles. 
























14 


C. MS-scheme 

In order to compare lattice values to experimental results one must convert to a universal renormalization 
scheme and use a reference renormalization scale. Typically one chooses the MS-scheme at a scale ^ of 2 GeV. 
For the conversion from the RI' to the MS scheme we use the intermediate Renormalization Group Invariant 
(RGI) scheme, which is scale independent and relates the RI' and MS results are follows: 




RGI 


c = AZ^^' in) = Z^^^(2 GeV) GeV). 


zRI'/ 


(50) 


The conversion factor can be read from the above relation: 


Z^^(2 GeV) = 2 GeV) (/i), 2 GeV) = 


AZ“S(2GeV) 


(51) 


where the scheme dependent quantity AZg(/x) can be expressed in terms of the /3-function and the anomalous 
dimension, yg = of the operator O (for definitions see Appendix El) : 




TO 

■ 200 


X ■'nw) 


7o 


(52) 


We employ the 3-loop approximation which simplifies to: 


AZ$i^l) = 2/3o 


IGtt^ 


2/3o 


1 -f 


g‘^(M)^ /3 i70 - Poll 
IGtt^ 2/3o 


(53) 


g (g) -20012 + /3o(7i (2/3i -h 7f) -h 2/327o) - 2/3o/3i7o(/3i -h 7f) -F 0 hl 


(1071^ 




where the coupling constant, is needed in both the MS and RI' schemes; their expressions coincide to 

three loops and read [3l| 


g 


MS.RI' 


(g) 


IGtt^ 


/3i log T 1 0f log^ L - 0f\ogL + 0200 - 0\ 


3-loop 00 L 0Q 0Q 


£3 


), T = log 


A2_' 

MS 


(54) 


For Aj^ we employ the value 315 MeV and 296 MeV for Nf=2 


and Nf=4: [^, respectively. 


VI. RESULTS 

In this section, we present our results for the renormalization functions in the MS scheme at a scale of 2 GeV. 
The final data correspond to the non-perturbative values after subtracting the lattice artifacts to 0{g^ a°°). 
Although, with the exception of .Zp, the dependence of the RFs on the pion mass is not significant, we nevertheless 
perform a chiral extrapolation of the RFs using data at the same 0 and Nf ensembles obtained for different pion 
masses as discussed in subsection IV Al 

As can be seen in Figs. [11]-[Ml the 0[g'^ a°“)-subtracted RFs (magenta diamond points) have a mild dependence 
on {ap0, which is removed by extrapolating to zero, using the Ansatz 

Zo(ap) = 4")+4')-(ap)2, (55) 

where Z^'^ corresponds to our final value on the renormalization functions. To extract the RFs reliably one needs 
to consider momenta in the range Aqcd < p < l/a. We relax the upper bound to be ~ 4/a to 7/a, which is 
justified by the weak dependence of our results on {ap)^. Therefore, for each value of 0 we consider momenta 
{apY > 2 for which perturbation theory is trustworthy and lattice artifacts are still small enough. 

From our analysis we find that the data for the RFs depend not only on (ap)^, but also on the directions of 
the momentum. Noting that, for democratic momenta (in all directions, not only spatial) the value of p4/p2^ 
equals 0.25, we find empirically that data produced on momenta with the ratio of Eq. (H^ . p4/p2^, being > 0.4, 


^ Sign differences in some terms of Eq. |[54l l compared to Ref. m are related to alternative definition of the /3-function 
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have a behavior that deviates from the general (ap)^ curve. The choice of such a momentum ratio as a criterion 
is justified by the fact that such Lorentz non-invariant contributions appear in the perturbative computation at 
higher orders in the lattice spacing (e.g., p4/p2 for 0{a^)). Thus, high values of this ratio are an indication of 
large lattice artifacts from higher loops. As an example, we demonstrate Zp^ in Fig. [15] including the data obtained 
at momenta satisfying P = p4/p2 > 0.4. These are shown by the filled blue circles and filled green diamonds 
corresponding to unsubtracted and a°“)-subtracted data, respectively. As can be seen, the filled symbols 
have different behavior than the open symbols. Although the subtraction of one-loop lattice artifacts reduces the 
difference, the higher order artifacts are not negligible. The data for these momenta have been excluded from 
the final analysis of all RFs. A similar study is presented in Refs. [3 [3. 
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FIG. 15: Za as a function of the renormalization scale. Filled blue circles and filled green diamonds correspond to 
unsubtracted and a°°)-subtracted data using momenta with P = pAjp'l > 0.4. 


While statistical errors are very small, a careful investigation of systematic errors is required. A small systematic 
effect comes from the asymmetry of our lattices, both because they are larger in their time extent and because of 
the antiperiodic boundary conditions in the time direction. To address this issue, we average over the different 
components corresponding to the same RFs, for instance Za is defined as: 

Zk = (^A + + ^a) (56) 

where the upper index indicates the Dirac matrix used as current insertion corresponds to insertion 75 7i). 
In addition, remaining systematics are automatically removed by the subtraction of the 0{g^ a°°) terms. The 
largest systematic error comes from the choice of the momentum range to use for the extrapolation to a?p^ = 0. 
One way to estimate this systematic error is to vary the lower or/and upper range used in the fit. Another 
approach is to fix a range and then eliminate a given momentum in the fit range and refit. The spread of the 
results about the mean gives an estimate of the systematic error. In the final results we give as systematic error 
the largest one from using these two procedures, which is the one obtained by modifying the fit range. 
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FIG. 16: Za (upper) and Zv (lower) as a function of the renormalization scale. The arguments a, b of Z{a,b) correspond 
to Nf and /3, respectively. From top to bottom, the data correspond to increasing the lattice spacing. Non-perturbative 
values are shown with the black circles, 0{g^ a^)-subtracted with the green crosses and 0{g^ a°“)-subtracted with magenta 
diamonds. 
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Fig. [THl corroborates that the magnitude of the lattice artifacts depends not only on the action parameters, but 
also on the operator under study, as can be seen for Zp^ and Zy shown for different values of the coupling constant. 
Since both Zy and Zp^ are scale independent, one expects a flat behavior as a function of the renormalization 
scale, {afj,)^ = (ap)^. However, the non-perturbative data before subtraction of the lattice artifacts is carried 
out, exhibit a non-zero slope, which becomes negligible once the 0[g^ a°°) terms are subtracted. For a proper 
comparison, we have kept the y-axis the same as the lattice spacing is increased. In the case of Za we find 
that the 0{g^ a?) terms computed for all iV/=4 gauge conhgurations, despite being the leading contributions, 
underestimate the total one-loop lattice artifacts, 0{g^ a^). Our analysis shows non-negligible lattice artifacts 
between 3 — 6% for momenta in the range [2,4]. Nevertheless, for the Nf=2 case, the total one-loop lattice 
artifacts are very small (0.1 — 2% for (ap)^ : [2,4]) which may be attributed to the inclusion of the clover term. 
One also observes that the 0{g‘^ a^) terms are no longer reliable, possibly due to the fact that they are polynomial 
functions of Csw (2 — 8% for the aforementioned momentum range). This fact is an evidence that the addition of 
the clover term in the twisted mass action suppresses lattice artifacts. This is also observed for other quantities 
besides renormalization functions, such as in the isospin splitting in the A-system [H, 

For Zv, on the other hand, we find that for all ensembles analyzed in this work, there are negligible one-loop 
artifacts beyond 0{g‘^ a?), as can be seen in the lower panel of Fig. [ini From our study we find that lattice 
artifacs are current-dependent and can be identified a posteriori, from the perturbative computation. This is also 
confirmed by examining the results using the Nf=2 ensembles with the clover term shown in Figs. [T71- 1231 where 
for Zp, the 0{g'^ a?) and 0{g^ a°°) are almost equivalent, especially for (ap)^ < 5; This is not the case for the 
other RFs shown in Figs.ITTl-l^ The Nf=2 are the most recent gauge configurations produced by ETMC, which 
are currently being used for hadron structure studies and thus the values of the RFs are needed to renormalize 
the hadron observables [13 . Figs. [171 -051 correspond to the RFs upon conversion to the MS scheme at 2 GeV 
and are plotted against the initial renormalization scale, (ap)^. 
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FIG. 17: Renormalization of the fermion field for Nf=2 twisted mass clover-improved fermions. The data correspond to 
the MS scheme at a reference scale of 2 GeV and are plotted against the initial renormalization scale, {ag)^ = (ap)^. 
Black circles (magenta diamonds, green crossed) denote the unsubtracted {0{g^ (-subtracted, 0{g^ a^)-subtracted) 
non-perturbative data. 
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FIG. 18: The renormalization function of the scalar operator. The notation is the same as that of Fig. 1171 
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FIG. 21: The renormalization function of the tensor operator. The notation is the same as that of Fig. 1171 
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FIG. 22: The renormalization functions of the one-derivative vector and axial operators. The notation is the same as that 
of Fig. \T7i 



FIG. 23: The renormalization functions of the one-derivative tensor operators. The notation is the same as that of 

Fig. El 

Comparing, for instance, Zy and Za computed using the Nf=A and Nf=2 ensembles (see Fig. [HI), we find 
that the 0{g'^ a°°) lattice artifacts for the Nf=2 ensembles are smaller and lead to good quality plateaus. The 
extrapolation to (ap)^ —0 is performed using the 0{g‘^ a°°)-subtracted data and for momenta with {apY > 2 
which is the range of interest. The data display a very small slope thus leading to a good determination of the 
continuum value. For Zp and Zp/Zs there is a stronger dependence on {apY up to (apY ~ 2 — 3. Thus, for Zp 
and Z-pjZ^ we use a different interval for the (apY 0 ht. For instance, in the Nf=2 ensembles plotted here, 
we ht in the range [4,7] for Zp and Zp/Zs, and in the range [2,7] for the remaining RFs. The systematic errors 
due to the choice of the ht, are computed by taking the difference in the values of Z^^ (see Eq. (l55l) l extracted 
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from these ranges and the range [3,5]. 

In Table Hill we give our final chiral extrapolated values of from a°“)-subtracted data (e.g. for Nf=2 
the hlled diamond point in Figs. [16]-[23]). Zp and Zp/Zs are obtained only at /3=2.10 (for both Nf=2,A) where 
we have enough ensembles for the pion-pole subtraction; the corresponding results are presented in Table IIVI 
Some of these RFs have been computed in Ref. [2l|; the small differences observed between our results and 
those of Ref. [2l| can be attributed to two factors: a) the different method for calculating non-perturbatively the 
vertex functions (see footnote 1) and b) the different analysis procedure since the authors use the subtraction 
of 0(ffboosted lattice artifacts, which are, in general, larger than the 0{g^ a°°) terms, leading, in some cases, 
to lower estimates. We have checked that if we apply the a^)-subtraction to our non-perturbative 

estimates, we are in agreement with the results of Ref. [2l| . Since the difference between our approach and that 
of Ref. is the treatment of lattice artifacts, both sets of results should agree in the continuum limit, a —>■ 0. 
This is indeed the case, as demonstrated in Fig.[24]for the vector and axial RFs. Our results are also compared to 
those of Ref. [S^. We find that upon taking the continuum limit, the differences between the two works become 
very small, or compatible with zero. We note that in Ref. [s^ the subtraction of lattice artifacts is performed via 
a hypercubic removal procedure. In addition, they use general momenta that are not restricted to democratic or 
near democtratic. 


RFs 

iV^=2,/3=2.10 

Ar^=4,d=2.10 

Ar^=4,d=1.95 

iV^=4,/3=1.90 

'7MS 

0.8366(2)(7) 

0.7822(4)(4) 

0.7835(2)(25) 

0.7480(6)(11) 

'7IVIS 

0.6606(9)(18) 

0.7143(9)(216) 

0.7342(1)(21) 

0.7835(2)(17) 

Zv 

0.7565(4)(19) 

0.6651(2)(5) 

0.6298(5)(29) 

0.6015(2)(4) 

Zk 

0.7910(4)(5) 

0.7744(7)(31) 

0.7556(5)(85) 

0.7474(6) (4) 

^MS 

Z/'p 

0.8551(2)(15) 

0.7875(9)(15) 

0.7483(6)(94) 

0.7154(6)(6) 

/7MS 

^DVl 

1.1251(27)(17) 

1.0991(29)(55) 

1.0624(108)(33) 

1.0268(26)(103) 

/7MS 

^DV2 

1.1396(16)(13) 

1.1398(37)(91) 

1.1209(61)(32) 

1.0676(44)(190) 

/7MS 

^DAl 

1.1494(9)(99) 

1.1741(42)(173) 

1.1255(27)(328) 

1.1151(51)(197) 

/7MS 

^DA2 

1.1357(20)(205) 

1.1819(47)(147) 

1.1555(36)(289) 

1.1170(54)(223) 

/7MS 

^DTl 

1.1377(160)(13) 

1.1562(32)(7) 

1.1218(106)(44) 

1.0777(37)(122) 

/7MS 

^DT2 

1.1472(121)(48) 

1.1822(59)(118) 

1.1727(121)(73) 

1.0965(90)(278) 


TABLE III: Our final values of the renormalization functions. The scheme and scale dependent RFs are given in MS at 2 
GeV. The number in the first parenthesis is the statistical error, while the number in the second parenthesis corresponds 
to the systematic error obtained by varying the fit range in the (ap)^ 0 extrapolation. 


RFs 

lV/=2,d=2.10 

Ary=4,/3=2.10 

/7MS 

Zp 

0.5012(75)(258) 

0.5468(15)(176) 

Zp [Z's, 

0.7016(141)(113) 

0.7036(23)(195) 


TABLE IV: Our final values for .Z|?®(2GeV) and ZpjZs- The number in the first parenthesis is the statistical error, 
while the number in the second parenthesis corresponds to the systematic error obtained by varying the fit range in the 
(op)^ —>■ 0 extrapolation. 
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FIG. 24: Difference of Zv and Za computed in this work and in Ref. (Method 1) for the A^/=4 case. The data are 
plotted against the lattice spacing, and their extrapolation to the continuum limit is shown with a dashed line. Solid 
points correspond to the limit a ^ 0; they are consistent with zero. 


VII. CONCLUSIONS 

We present results on the renormalization functions of the fermion field and fermion bilinears with up to one 
covariant derivative. The computation is performed non-perturbatively on several ensembles of Nf=A twisted 
mass fermions, as well as Nf=2 twisted mass fermions including a clover term. This work is a continuation of 
our renormalization program first addressed in Refs. Besides the analysis of the IV/=2 twisted mass 

clover-improved ensembles we have improved the procedure for the subtraction of lattice artifacts. The procedure 
that we adopt here for the perturbative computation of lattice artifacts is improved by taking into account not 
only leading order lattice artifacts, 0{g‘^ a^), but also contributions to all orders in the lattice spacing, 0{g^ a°°). 

The non-perturbative computation uses a momentum-dependent source and the RFs are extracted for all the 
relevant operators simultaneously. This leads to a very accurate evaluation of the RFs using only a small ensemble 
of gauge configurations (0(10)). The precision of the results allows us to reliably investigate the quark mass 
dependence, which is found to be very weak with the exception of Zp. Nevertheless, a linear extrapolation with 
respect to the pion mass squared is carried out in order to reach the chiral limit. For the renormalization function 
of the pseudoscalar operator, Zp, we find a quark mass dependence due to the pion pole, which we eliminate 
using a refined procedure that avoids introduction of artificially large errors. The procedure entails a suitable 
fit (Eq. (H5]) l to identify Zp directly from the constant term, ap, instead of subtracting the pion-pole term (see 
Eq. (HTlll. 

Our accurate non-perturbative results show that, although the lattice spacings considered in this work are 
smaller than 0.1 fm, lattice artifacts are not negligible in most cases, and are significantly larger than statistical 
errors. Thus, the subtraction of the 0{g^ a°°) perturbative contributions appear to improve significantly the 
determination of the RFs, by leading to a milder dependence of the RFs on {apY. Residual 0{a?'p^) effects are 
removed by extrapolating our results to (ap)^ = 0. For the scheme and scale dependent RFs, we convert our 
values to the MS scheme at a scale of 2 GeV. The statistical errors are, in general, smaller than the systematic 
ones. The latter are estimated by changing the window of values of the momentum used to extrapolate to 
a?p^ = 0. Our final values are given in Tables Hill - IIVI 
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Appendix A: /^—function and anomalous dimensions 

For completeness we provide in this Appendix the definition of the /3—function and the anomalous dimension 
of the operators studied in this work, which include up to one covariant derivative. To simplify the expressions 
we give the perturbative coefficients in the Landau gauge and in SU{3). 
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The perturbative expansion of the anomalous dimension in a renormalization scheme S is given as follows: 


o d , ^ 9 


IGtt^ 
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IGtt^ 


Similarily, the /3—function is defined as: 
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/3 (g) Po ^1(16^2)2 
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For the conversion from the RF to the MS scheme we use the three-loop expressions, to which the coefficients of 
the /3—function coincide and are given by [H, H^: 


Po 

Pi 

P2 


102 -flV/, 

2857 5033 
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Ar 325 

Nf H- N^f. 
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(A3) 

(A4) 

(A5) 


Below we give all necessary expressions to convert to the MS scheme, as well as the references from which they 
were taken (see also references therein). Dome signs and multiplicative numerical factors have been adjusted to 
match the definition of Eq. (EH. An upper index appears for scheme-dependent quantities, in order to denote 
the scheme that they correspond to. 

Quark field [4l| : 


Tensor l40l l4 
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(C 3 = 1.20206...) 

Scalar/pseudoscalar [4 
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One-derivative vector/axial (^. 1^: 
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One-derivative tensor 14511: 
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Appendix B: Application of the subtraction to 0{g^ a°°) in other ensembles 

As discussed in the main part of the paper in our previous works of Refs. [Hill we have applied a procedure of 
subtracting the lattice artifacts of 0{g^ a^). The values of the RFs are used to renormalize hadron quantities such 
as the axial charge and the quark momentum fraction, in order to compare them with other lattice discretizations, 
as well as with experimental data. For a fair comparison between renormalized matrix elements of the ensembles 
presented in this work and the ones given in Refs. we have updated the RFs of the latter publications by 

applying the subtraction procedure to one-loop and all orders in the lattice spacing, 0{g^ a°°). These correspond 
to tree-level Symanzik improved gauge action and Nf=2 twisted mass fermions at three values of the coupling 
constant corresponding to /3=3.90, 4.05, 4.20. Since the gluon action is different from the ensembles of Table HI 
and since employed momentum values are also different, a perturbative computation of the 0{g^ a°°) contributions 
was required on each ensemble in order to match its parameters, such as the coupling constant, the lattice size 
and the values of the renormalization scales. The new data on the Rfs are given in Table fVl 


RFs 

p = 3.9 

P = 4.05 

P = 4.20 

/ 7 MS 

0.769(1)(2) 

0.787(1) (3) 

0.783(1)(2) 

/ 7 MS 

0.791(2)(41) 

0.748(2)(31) 

0.754(1)(16) 

/ 7 MS 

Zp 

0.527(6)(70) 

0.517(2)(33) 

0.546(5) (33) 

Zp/Zs 

0.672(7)(60) 

0.700(3)(14) 

0.731(5)(25) 

Zv 

0.646(2)(2) 

0.681(2)(6) 

0.701(1)(4) 

Za 

0.769(2)(1) 

0.787(1)(1) 

0.791(1)(1) 

' 7 MS 

Zj'y 

0.758(2)(4) 

0.796(1)(3) 

0.814(1)(3) 

^MS 

^DVl 

1.028(2)(6) 

1.080(2)(11) 

1.087(3)(12) 

/ 7 MS 

•^DV2 

1.064(4)(4) 

1.123(4)(10) 

1.130(4)(4) 

' 7 MS 

•^DAl 

1.106(3)(8) 

1.157(4)(10) 

1.150(4)(15) 

/ 7 MS 

^DA2 

1.102(5)(7) 

1.161(4)(13) 

1.164(3)(6) 


TABLE V: Updated results on the RFs of Refs. [l^ . [l^ {Nf=2, /3=3.90, 4.05, 4.20, Csw = 0) using the subtraction 
procedure to 0{g^ a°°). 
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